From extended phase space dynamics to fluid theory 
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We derive a fluid theory for spin- 1/2 particles starting from an extended kinetic model based on a spin-projected density 
matrix formalism. The evolution equation for the spin density is found to contain a pressure-like term. We give an 
example where this term is important by looking at a linear mode previously found in a spin kinetic model. 
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I. INTRODUCTION 

Quantum plasma physics is an active field of research and 
the methods have been applied to investigate a wide range 
of different effects. Some recent examples are quantum ion- 
accoustic waves^ Jeans instabilities in quantum plasmas^, 
trapping effects^, magnetization by photons'* and relativistic 
effects^'^. Typically, the quantum hydrodynamic equations 
are derived starting from the Schrodinger equation and mak- 
ing the Madelung ansatz^. However, a method that more 
closely resembles the classical case is to use kinetic equations 
as a starting point, see Ref. [7] for a comparison between the 
different approaches. The field of quantum kinetic theory** 
in many ways started with the ambitions of Wigner, as pre- 
sented in Ref. [9], to bridge the gap between classical Liou- 
ville theory and statistical quantum dynamicsi^— . Thus, the 
development of quantum kinetic theory was partly due to an 
interest in obtaining a better understanding of the quantum- 
to-classical transition'^. However, another important aspect 
of quantum kinetic theory is as a computational tool for, e.g., 
quantum plasmas'"*, condensed matter systems'^''-, and, in 
general, quantum systems out of equilibrium'^, and in that re- 
spect shares many commonalities with quantum optics'**. 

Systems where the electron spin plays an important role is 
a growing field of research. One reason for this is the promise 
of novel electronic devices. Furthermore, the utilization of 
spin may also play a key role in the yet-to-come implementa- 
tion of quantum computin g '^•^'^ . There is also a growing in- 
terest in quantum effects in laser generated plasmas^'"^^ and 
other high energy density environments^ Thus, for such sys- 
tem, more encompassing modeling is required as compared 
to the condensed matter system used in modeling of, e.g., 
electronic devices. Spin effects has also been studied in the 
context of thermonuclear fusion plasmas by Cowley, Kulsrud 
and Valec^^, who investigated spin depolarization rates. Their 
idea was to start from a kinetic theory for spin- 1/2 particles 
and calculate moments of this to obtain macroscopic quanti- 
ties. The distribution function in their theory was defined on 
an extended phase-space where the spin degree of freedom 
was described by vector The kinetic equation, however, was 
not derived from a fully quantum mechanical theory and this is 
a point that later was criticized by Zhang and Balescu^^, who 
derived a fully quantum mechanical equation for spin- 1/2 par- 
ticles. Zhang and Balescu also argued that the spin cannot be 
described by an extended phase-space and thus their kinetic 



approach is based on a Wigner matrix. However, it is possible 
to describe the spin through an extended phase-space, see^- 
for a review of quasi-distribution function for spin. In an sub- 
sequent paper^** Balescu and Zhang then went on to derive a 
fluid moment theory based on their kinetic theory. This the- 
ory was derived to leading order in Ti which means that the 
dynamics due to the spin dipole force is not retained. 

Recently, we have derived a scalar kinetic theory for spin- 
1/2 particles^^'^'^, where the spin is described by a unit vec- 
tor in the distribution function. This is natural extension of 
regular quantum kinetic theories^' as well as fluid models 
with spin (see e.g. ll32ll and lHH). This theory has the ad- 
vantage of being of close mathematical resemblance to the 
Vlasov theory for a classical system, where the spin degree of 
freedom is seen to correspond to a classical magnetic dipole 
moment. This may guide the intuition when interpreting the 
corresponding fluid moment theory. It may also make classi- 
cal numerical codes amendable and adaptable for quantum ki- 
netic simulations. The scalar quantum kinetic theory^" is seen 
correspond to the kinetic equation of Cowley, Kulsrud and Va- 
leo (in the collisionless limit) to leading order in h. Our start- 
ing point in this paper is a kinetic equation valid to order /i^, 
which will then keep not only precession but also the dipole 
force acting on the spin. 

Spin fluid models have also been applied in other contexts, 
for example considering ultra cold trapped atomic gases^"*. In 
this case, the atoms interact via two hyperfine states and these 
can be viewed as the two spin states, hence the important fea- 
ture of such a theory is the spin-spin interaction. In the theory 
considered here, however, we will not deal with two-particle 
spin-spin interactions, but rather use a mean-field approach. 

In the presence of a vector potential, the gauge issue 
makes itself reminded, since a gauge transformation A^(x) ^> 
+ dfiXi^) of the four- vector potential gives the trans- 
formation v/(jc) '^{x)&K^[iex{x) /he] of the wave function. 
Thus, in order to obtain a gauge invariant quasi-distribution, 
the ansatz for the Wigner function has to be modified. There 
is infinitely many ways to do this and external criteria have to 
be applied in order to choose a specific modification^^. Thus, 
there are different 'standard routes' for the modification due 
to the vector potential, e.g. using a Wilson looj>^. In quantum 
kinetic theory the standard way is, to the best of our knowl- 
edge, best represented by Stratonovich and others^^-^**, who 
choose to take the gauge modification of the wave function 
into account by adding a line integral of the vector potential 
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to the canonical momenta in the definition of the Wigner func- 
tion. However, the corrections due to the gauge choice enters 
the phase space dynamics at higher order in the Planck's con- 
stant h. As we in this paper will be mainly interested in the 
lowest order spin corrections to the fluid equations, we intend 
not to make use of the full gauge operators that comes from 
short length scale (compared to the characteristic de Broglie 
wavelenght) contributions. The fluid equations based on the 
full kinetic theory has been investigated recently in the spin- 
less case^2ilS. 

In Section 2 we derive the fluid moments of the spin kinetic 
equation and discuss the derived model. The equations are 
then coupled via Maxwell's equations and we consider linear 
modes of the system in Section 3. We also briefly discuss the 
closure problem. In Section 4 we draw our conclusions and 
discuss some future applications. 



II. MOMENTS OF THE SPIN KINETIC MODEL 

In a recent paper— we derived a kinetic model for a par- 
ticle with spin- 1/2. In this theory the system is described 
by a quasi-distribution function /(x,v,s,f) which is a com- 
bination of the Wigner function^ in phase-space variables x 
and V and a transformation for the spin variable, the unit 
vector s. It has the property that the marginal distribution 
n(x, s) = / d^vf yields the probability density to find a particle 
in position x with spin in the direction given by s. Similarly 
n(v,s) = / d^xf gives the probability to find a particle with a 
velocity v and spin in the s-direction. A third important prop- 
erty is that the magnetization M due to the spin is obtained 
by 



M(x,f) = 3jti dQ.sf{x,\,s,t), 



(1) 



where dQ. — d^vd^s and jx is the magnetic moment of the 
particles. Here we will consider electrons and hence ji = 
— {g/ 2) eh /{2m) where e is the elementary charge, m is the 
electron mass and h is Planck's constant divided by 2n. The 
g-factor for an electron is modified from the value g — 2, pre- 
dicted by the Dirac theory, to § « 2.0023 when one takes vac- 
uum fluctuations into account, see Ref.*^ for a recent mea- 
surement. Note the factor 3 which is related to the non- 
commutation of the different spin components. 

If we consider scale lengths longer than the characteristic 
de Broglie wavelength the evolution equation for the quasi- 
distribution function is found to bo^ 

5r/ + v-V,/+^(E + vxB)-V„/ 
m 

+ ^ V, [(s + V,) . B] • V„/+ ^ (s X B) • V,/ = 0. 
m n 



(2) 

The terms proportional to /i are terms that arises due to the 
spin, the first term is the dipole force for a particle with mag- 
netic moment /is, which has been modified due to the quan- 
tum mechanical nature of the spin, by the term (/i/m)Vv(B • 



Vj) ■ Vv/. The last term on the second row is due to spin pre- 
cession. To obtain a macroscopic theory for the system we 
define the fluid moments 



d£if, 



Ui = - I davif, 



P. J - 
S,= 



mj da{vi-Ui){vj-Uj)f, 



ILij - m / dD.{-isi - Si){vj - Uj)f, 



da{3si~Si)ivj-Uj){v,-U,)f. 



(3) 
(4) 
(5) 
(6) 
(7) 
(8) 



Note that these definitions are in line with the results of Ref. 
ll28ll . However, our definition of the mixed spin-velocity mo- 
ment O contains only the thermal contribution. The first three 
moments, Eqs. (O to (|5]l, are the density, velocity and pres- 
sure. The fourth moment, Eq. yields the spin density per 
volume. The two following moments are mixed spin and ve- 
locity moments. They have no straightforward physical inter- 
pretation, but we will see that the spin-velocity moment (|7| 
will act as a pressure-like term in the evolution equation for 
the spin density. 

Using the definition for the different moments above it is 
straightforward to derive evolution equations using Eq. (|2|i. 
The evolution of the density moment is given by the continuity 
equation 



The fluid velocity is evolving according to 
DUi q 



Dt 



m ^ dxi 



(9) 



(10) 



dxj 



where summation over repeated indices is implied and the 
convective derivative D /Dt = d /dt + Ujd /dxj has been intro- 
duced. The second term on the right hand side is a magnetic 
dipole force due to the spin. The pressure moment is evolving 
according to 



Dt 



'^ik jk '^ij 

dxji dxji axit 

q q 

m m 



dUi 



H„ dBt U dBk 

^ik^ \ ^ik^ 

m dXj m dXi 



dQ, 



dxv 



and as expected has some extra terms due to the spin (cf. 
The evolution equation for the spin density is given by 



DSi 



1 dZi 



ZZl = ^£.. S B ^ "^'^ 

Dt h '■''^ ■' ^ nm dxi 



(11) 
39i). 

(12) 



Comparing this with Eq. ( fTOl i we see that this mixed spin- 
velocity tensor is acting as a pressure term on the spin mo- 
ment. Exactly which physics is captured by this term is mainly 
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a subject of further research. In this paper we give an exam- 
ple where this term is important. Finally, the spin-velocity 
moment evolution equation is 



OS, 



Dt 



!1 = _E,.^_E,,M_P.,^ 
'■' dxj^ dxji dxk 

q „ 2/i 

H Ejkl^ikBl + —£ikl2lkjBl 

m n 

SB, dBk dAijk 

+lln- linSiSk^ 

axj axj axj^ 



(13) 



It is instructive to compare these fluid equations with the ones 
obtained in Ref. |32] which is based on a Madelung ansatz 
of the wave function. The equations obtained there contain 
additional terms which are explicitly quantum mechanical al- 
ready in the velocity and spin moment equations, for example, 
the Bohm potential and various spin-spin correlation terms. In 
the current paper these terms are missing for two different rea- 
sons. Some of the terms are 'hidden' in higher order moments 
(see lHH), and others are missing since our starting point here 
was the long scale length limit of the kinetic equation. 



III. LINEAR MOTION PERPENDICULAR TO THE 
MAGNETIC FIELD 

We now consider an electron-ion plasma which is subject 
to an external magnetic field in the z-direction, Bq = BqI. If 
we consider sufficiently fast dynamics the ions can be taken 
as a stationary neutralizing background with number density 
«o. To obtain a closed set of equations the closure conditions 
Qijk — and Aiji, — are assumed, appropriate for low tem- 
peratures. For perturbations of the systemn = «o + «i, v = Vi, 
Pij = S = So + Si, E = El and B = Bo + Bi, with 

small deviance from the equilibrium (denoted by subscript 0) 
we may linearize the equation of motion for the electrons Eq. 

dinii-dni. 

We will assume a thermal background distribution /q. 
The appropriate form of this distribution function is (for 
sufficiently large chemical potential) a generalization of a 
Maxwell-Boltzmann distribution given by^ 



no 



(27r)V2v: 



-vV(2v|) 



l+tanhf^Vose. 



(14) 



Using this distribution the equilibrium pressure moment is 
Po ~ nokBT and the equilibrium polarization density is 



(15) 



A. Two-fluid model 

The Eqs. (fT0]l-(fT3Tl constitutes spin-fluid equations for a 
single electron-fluid. However, since the kinetic equation is 
linear in /, it is straightforward to divide the electron fluid 
into two parts depending on the initial spin state, e.g. up- or 
down relative to the magnetic field. The Eqs. (fT0li-(fT3]l will 
then apply to two sets of electron fluids, which differ by their 
initial value of the spin moment (|6]l. Such a division has been 
succesful in deriving the rate of spin-polarization due to the 
spin-dependent part of the ponderomotive forced. In case 
simplified spin-fluid models without the velocity-spin tensor 
are used (e.g. Refs. |33] and [41]), such a division is neces- 
sary, in order to capture the basic physics of how the magnetic 
dipole force can induce a spin-polarization in a plasma'*^. If 
the model is more refined and E,y is kept, on the other hand, a 
single fluid model will suffice to cover this effect. Therefore 
the choice of a one-fluid or two-fluid model is to a large ex- 
tent a matter of taste, provided E,y is kept. The advantage with 
the one-fluid model is that it reduces the number of dependent 
variables, since only a single electron-fluid is considered. The 
advantage with a two-fluid model (with division of the elec- 
tron fluid made with respect to the initial spin state) is that 
much of the physics can be captured in the simpler moments 
(|4]l and (|6]l, that are easier to interprete physically, rather than 
in the more complicated tensor (|7|l- It should be stressed that 
more highly nonlinear and/or thermal effects may require an 
evolution equation for even higher order moments, i.e. for 
Aijii and higher. 



where the last approximation is valid when fiBo <C ksT. Note 
that for electrons /i < so that the spins tend to align anti- 
parallel to the magnetic field. 

To find linear solutions we Fourier decompose with wave 
vector k — kx. Looking for solutions with the electric field 
polarization along the z-axis E = E-z we can in the end ver- 
ify that Gxz and a,,- are zero and hence it is only necessary to 
keep the av- component of the conductivity tensor For con- 
venience we here use an effective conductivity tensor where 
the total current j,, including both the free and the magnetiza- 
tion current, is related to the electric field by j,- = OijEj. We 
obtain, after some algebra, the dispersion relation 



(0 



'k^c^^co^- 



IrnksT 



(Or 



CO - COcg 

CO - COcg co + co, 



COcf, - COc 
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Co} I6m^ 

COc + (Ocg 
CO - COcg - Oc 

(Peg - COc 

CO + COcg -(Oc (0 + (Oca 



CO + COcg 

2 COcg 



■eg 



CO 



COc 



COcg + (Oc 
COcg 



COc 



(16) 



where we have kept only the terms that survive in the limit 
r — > 0. As can be seen from this there is a resonance at 
CO — COcg, however, we know from the corresponding kinetic 
theory— that we also have cyclotron modes with resonace at 
CO = COc- So, unless the frequency resolution is very high it is 
difficult to separate this spin resonance from the classical res- 
onance. Instead we consider frequencies sufficiently close to 
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the differences co sa A© ~ a)c{g/2 — l), the dispersion relation 
is then 



co^ - k^c^ = col - 



cot 



Aco 



co} 16ot2 CO — Aco 



(17) 



This dispersion relation is the same as obtained from the full 
kinetic theory in ll29tl . If one is interested in thermal effects it 
is necessary to keep the terms from the third order mixed mo- 
ment A,-^^, since thermal terms will be introduced in the second 
order moment which are cancelled by a contribution from the 
next higher order This seems to be a general complication 
with quantum fluid theory in the linear limit^^ ''^''. Since the 
linear solutions are obtainable directly from the kinetic theory, 
perhaps the main strength of using a fluid theory is when con- 
sidering non-linear motion, where the kinetic theory quickly 
becomes unmanageable. 



IV. SUMMARY AND DISCUSSION 

We have in this paper derived the fluid equations for spin- 
1/2 particles. It is found that important effects occur due to 
the correlation between spin and velocity, as covered by the 
velocity-spin tensor E,^. The closure problem of the hierarchy 
of fluid equations has been pointed out. We note that it is hard 
to know a priori how to close the fluid system in a given situ- 
ation. For a low-temperature linear problem, we have shown 
that it suffice to drop moments Ajji^ and higher. However, the 
real need for a fluid model occurs in inhomogeneous and/or 
nonlinear problems, when kinetic theory becomes increas- 
ingly cumbersome to apply both analytically and numerically. 
To what extent dropping higher order moments is adequate in 
such scenarios is an open question. An important area of fu- 
ture research is thus to investigate the applicability of the fluid 
model to more complicated problems. 
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